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Barber's second self creation theory with bulk viscous fluid source for an LRS Bianchi type-I metric 
is considered by using deceleration parameter to be constant where the metric potentials are taken as 
function of x and t. The coefficient of bulk viscosity is assumed to be a power function of the mass 
density. Some physical and geometrical features of the models are discussed. 

1. Introduction 

Several modification of Einstein's general relativity have been proposed and extensively 
studied so far by many cosmologists to unify gravitation and many other effects in the 
universe. Barber has produced two continuous self-creation theories by modifying the 
Brans-Dicke theory and general relativity. The modified theories create the universe out 
of self-contained gravitational and matter fields. Brans-Dicke a theory develops Mach's 
principle in a relativistic framework by assuming interaction of inertial masses of funda- 
mental particles with some cosmic scalar field coupled with the large scale distribution of 
matter in motion. However, Barber 111 has included continuous creation of matter in these 
theories. The universe is seen to be created out of self-contained gravitational, scalar and 
matter fields. However, the solution of the one-body problems reveal unsatisfactory charac- 
teristics of the first theory and in particular the principle of equivalence is severely violated. 
The second theory retains the attractive features of the first theory and overcomes previous 
objections. These modified theories create the universe out of self-contained gravitational 
and matter fields. Recently, Brans u has also pointed out that Barber's first theory is in 
disagreement with experiment as well as inconsistent, in general, since the equivalence 
principle is violated. Barber's second theory is a modification of general relativity to a 
variable G-theory. The second is an adaptation of general relativity to include continuous 
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creation and is within the observational ambit. In this theory the scalar fields do not directly 
gravitate, but simply divide the matter tensor, acting as a reciprocal gravitational constant. 
It is postulated that this scalar field couples with the trade of energy momentum tensor. In 
view of the consistency of Barber's second theory of gravitation, we intend to investigate 
some of the aspects of this theory in this paper. 

Several cosmologists have studied various aspects of Robertson-Walker model in Bar- 
ber's second self-creation cosmology with .perfect fluid satisfying the equation of state 
p = (7 — l)p, where 1 < 7 < 2. Pimentela and SolengB have discussed the Robertson- 
Walker solutions in Barber's second self-creation theory of gravitation by using power law 
relation between the expansion factor of the universe and the scalar field. Carvalho □ stud- 
ied a homogeneous and isotropic model of the early universe in which parameter gamma of 
'gamma law' equation of state varies continuously with cosmological time and presented 
a unified description of early universe for inflationary period and radiation-dominating era. 
Singh, I Reddy, Bi and Reddy et al have presented Bianchi type space-times solutions 
in Barber's second theory of gravitation. Reddy and Venkateswarlu present Bianchi type 
— VIq cosmological solutions in Barber's second theory of gravitation both, in vacuum as 
well^s in the presence of perfect fluid with pressure equal to energy density. Shanthi and 
Rao t3 studied Bianchi type II and III space-times in second theory of gravitation, both in 
vacuum as well as in presence of stiff-fluid. Recently, Shri Ram and Singh 113 have dis- 
cussed spatially homogeneous and isotropic R-W model of the universe in Barber's second 
self-creation theory of gravitation in the presence of perfect fluid by using 'gamma-law' 
equation of state. Barber's second self creation theory with perfect fluid source for an LRS 
Bianchi type-I metric is considered recently by Pradhan and Vishwakarma li-3 by using de- 
celeration parameter to be constant where the metric potentials are taken as function of x 
and t. 

Most studies in cosmology involve a perfect fluid. However, observed phenomena 
such as the large entropy per baryon and the remarkable degree of isotropy of the cosmic 
microwave background radiation, suggest that we should analyse dissipative effects in cos- 
mology. Further, there are several processes which are expected to give rise to viscous 
effects. These are the decoupling of neutrinos during the radiation era and the decoupling 
of matter and radiation during the recombination era. Bulk viscosity is associated with the 
GUT phase transition and string creation. It is known that the introduction of bulk viscosity 
canavpid the big bang singularity. For further discussions, see the references in Pradhan et 
al.0-@ 

For simplification and description of the large scale behaviour of the actual universe, lo- 
cally rotationally symmetric [henceforth referred as LRS] Bianchi I spacetime have widely 
studied. When the Bianchi I spacetime expands equally in two spatial directions it is 
called locally rotationally symmetric. These kinds of models are interesting because Lidsey 
@ showed that they are equivalent to a flat (FRW) universe with a self-interacting scalar 
field and a free massless scalar field, but produced no explicit example. Some explicit solu- 
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tions were pointed out in Refs. EH Motivated by these above arguments, in this paper, we 
have investigated spatially flat and non-flat bulk viscous LRS Bianchi type-I cosmological 
models in Barber's second theory of gravitation. 



2. Field Equations 

We consider an LRS Bianchi type-I spacetime 

ds 2 = dt 2 - A 2 dx 2 - B 2 (dy 2 + dz 2 ) (1) 

where A = A(x,t),B = B(x,t). The field equations in Barber's second self-creation 
theory are 

Rij - 2 9ijR = 0^ ^ 

and 

Ucj) = 4ttAT (3) 

where □</> = </)', is the invariant d' Alemberian and T is the trace of the energy momentum 
tensor. A is coupling constant to be determined from the experiment (|A| < 0.1). In the 
limit A — * this theory approaches the standard general relativity theory in every respect 
and G = 4 . The stress-energy tensor in the presence of bulk stress has the form 

Tij = (p+ p)UiUj - pg i0 (4) 

where 

p = p — and u l Ui (5) 

Here p, p, £ and u are, respectively, the energy density, isotropic pressure, bulk viscous 
coefficient and four- velocity vector of the distribution. Corresponding to metric (1), the 
four velocity vector Ui satisfies the equation 

'I, j = 1 (6) 
The Bianchi identities in contravariant form applied to Eq. (Q) are 

U ; , 7 '■' + WT? = (7) 

where W = Sncfr^ 1 and in general relativity W = — 8irG. A comma and a semicolon de- 
note ordinary and covariant differentiation, respectively. In a comoving coordinate system, 
the surviving components of the field Eqs. for metric (1) are 

2B B 2 B' 2 , , 

-B + W - A^ = P (8) 
.. B'A 

B — = (9) 
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Aqb 2Bql A'qb' 2B'4>' 



57T(p p 



(10) 
(11) 
(12) 



A B A 3 A 2 B 
Here and in following expressions a prime and a dot indicate partial differentiation with 
respect to x and t respectively. 

3. Solution of the Field Equations and Discussion 

In this section, we review the solutions obtained by Pradhan and Vishwakarma. 
after integration, yields 

& 



A = T 

where I is an arbitrary function of x 

Eqs. (§h and (IhJ), with the use of Eq. (115), reduce to 



(13) 



B d B 



B d (B 



B< dx\B 



B'dtKB 



I 2 



B 2 



BV 



B'l 



(14) 



If we assume as a function of x alone, then A and B are separable in x and t, Eq. ( p"4| ) 
gives after integration 

B = lS(t) (15) 
where S(t) is an arbitrary function of t. With the help of Eq. (|l5]), Eq. ( |l3| ) becomes 



A -J 7 s 

Now the metric ( 1 ) takes the form 

ds 2 = dt 2 - S 2 (t)[dX 2 + e 2X (dy 2 + dz 2 )} 



where X = lni. With the use of Eqs. (|15|) and (|16|), Eqs. (g), (|llj) and (]12j) yield 

25 S 2 1 o .-i 

~S + £2 ~ = ~ 8n ^ P 

3 35 2 n , , 
52 "55- = 8^ P 

7 305 30'? Z 2 d 



I'S 2 I'S 2 dx \ I 



8ttA 



(p - m 



For the sake of simplicity, if we assume <j> to be a function of £ only, then Eq. (£ 
use of Eqs. ( p"8| ) and ( p"9| ) reduces to 



30S* 2AS 



05 



5 



= 



(16) 
(17) 

(18) 

(19) 

(20) 
with the 

(21) 
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The function S(t) remains undetermined. To obtain its explicit dependence on t, one may 
have to introduce additional assumptions. In the following, we assume the deceleration 
parameter to be constant to achieve this objective i.e. 



55 
Is 2 



H + H 2 
H 2 



^(constant), 



(22) 



where, H — 5/5 is the Hubble parameter. The above equation may be rewritten as 

5 ,5 2 
5 + ^ = ° 

On integration Eq. (|23|) gives the exact solution 



S(t) 



[a(t - t )j < 1+b > when b ^ -1 



m 2 t 



when b = — 1 



(23) 



(24) 



where a, mi and TO2 are constants of integration and the constant to means the freedom of 
choosing the time origin. Using Eq. @ in Eqs. ©, @ and (fn|) lead to 



^0-^=^ + (26 -1)J? 2 



Jtt^p = 3(^j - ff 2 ) 



2\bH 2 = 



(25) 
(26) 
(27) 



3.1. Non-flat Models 

Case (i): 6 7^ —1. For singular models since 5(0) = 0, Eq. (^4j) leads to 

5 = mt^+^ 
Using Eq. (gf) in Eqs. ©, © and © yield 

= m-ff^ra LH'+M + C2e f (i-VI=isd' 



where 



Btt(p - = 
87rp = 3</> 



J (26-1) 

^(TTb) (l + 6) 2 i 2 



1 



1 



ki = 



mHTT+b) (1 + b) z t 2 

[26(4A + 3) - 3] 
4(1 + b) 2 



(28) 

(29) 
(30) 

(31) 
(32) 
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and ci and c 2 are arbitrary constants of integration. Therefore, the geometry of the universe, 
in this case, is described by the line-element 



ds 2 = dt 2 - m z tT^) [dX 2 + e 2X (dy 2 + dz 



(33) 



Here the Barber's scalar function qb is given by Eq. ( f29| ) and the corresponding physical 

parameters p and p are given by the Eqs. ( ^0| ) and ([H]) respectively. 

Thus, for given ^(i), we can solve the pressure p. It is standard to assume the follow- 
ing ad-hoc law 

m = ZoP n (34) 

If n = 1, Eq. ( [34] ) may correspond to a radiative fluid, whereas n_= | may correspond to 
a string dominated universe. However, more realistic models are based on n lying in 
the region < n < h. 
Model I (£ = £0) 

In this case, we assume n = 0. Further from Eq. (p4]), we obtain £ = £0 = constant and 
hence Eq.(|30[) becomes 



1 



(26-1) 



2£o 



8nt 2l m 2 t -rTh) (l + b) 2i (l + b)t 



(35) 



Model II (£ = £ p) 



In this case, we assume 11 = 1. Further from Eq. (|34j), we obtain £ = ^oP an d hence Eq. 
(DO) reduces to 



1 



P 



(2b -1) 

87rf 2l TO 2;(TTbT (l + b) 2 



9£ 



-{ 



1 



1 



(1 + 6)* (! + & ) : 



■}] 



(36) 



From Eq. (|Tj), we observe that the energy density pit) is decreasing function of time. As 
t tends to infinity, energy density will vanish. From Eq. ( p9| ) , it is observed that Barber's 
scalar function qb decreases as time increases and will vanish when t turns to infinity. The 
energy conditions E3 'E3 are satisfied provided m 2 < 1 and 6 = with positive c% and c 2 . 
Physical behaviour of the models: In case of a non-flat model when 6^—1, the Ricci 



scalar becomes 



R 



(1-6)* 



(37) 



m 2£(T+6) (1 + 6) 

It is observed from Eq. ([37]) that when * — > 0; (i) R — > 00 if 6 = 0, and (ii) i? — > 00 if 
6 > 1. The Eq. ( [37| ) also suggests that when * — ► ex); i? — > if 6 > 0. The scalars of the 
expansion and shear are given by 
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The model has singularity at t = 0. At t — ► 00, the expansion ceases. Here, ? = 0, which 
confirms the isotropic nature of the space-time which we have already obtained in Eq. (B3b. 



Case (ii): b = — 1. In this case using Eq. (E4J) in Eqs. (27), (E5j) and (E6J) give 



mi 



8tt(p + 3m 2 £) = 
87175 = 3</> 



3m 2 " 



mi 2g2m 2 * 



1 



mi 2 e 2m2t 



m,2 



where 



m 2 



2A +7 
4 



(39) 
(40) 
(41) 

(42) 



and C3 and C4 are arbitrary constants of integration. Therefore, the geometry of the universe, 
in this case, is described by the line-element 



ds 2 = dt 2 - m 1 2 e 2m * t [dX 2 + e 2X (dy 2 + dz 2 )] 



(43) 



In this case, the Barber's scalar function is given by the Eq. ([39|) and the corresponding 
physical parameters p and p are given by Eqs. (^0|) and (^lj) respectively. For ni2 = 0, the 
model reduces to a static radiating model with constant density and pressure. 
Model I (£ = £ ) 



In this case, we assume n = 0. Further from Eq. (34), we obtain £ = £o= constant and 
hence Eq. (^o|) reduces to 



1 



P : 



87T m 2 e 



2 O 2mof 



37713] - 3m 2 £o 



(44) 



Model II {£ = £ p) 



In this case, we assume n = 1. Further from Eq. (|34|), we obtain £ = £o/° an d hence Eq. 
( j40| ) reduces to 



- 3m 2 , - 9m 2 £o{ 



87r L mfe 2m2t 

Physical behaviour of the model: The Ricci scalar R is 



;}] 



(45) 



R = 2m 2 2 - 



1 



j 1 2 e 2m 2 t 



(46) 



It is easily observed from Eq. ( p^| ) that (i) when t — ► 0, i? — > (2m 2 2 — ^pr)> an d (ii) 
when i — > 00, i? — > 2m 2 2 . The expansion and shear scalars are 



6 = 3m 2 , (7 = 



(47) 
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The model represents an uniform expansion as can be seen from Eq. The flow of the 
fluid is geodetic as the acceleration vector fi = (0, 0, 0, 0). 



3.2. Flat Models 

Motivated by the results of the BOOMERANG experiment on Cosmic Microwave 
Background Radiation, l3 we wish to study a spatially flat bulk viscous cosmological mod- 
els. 

The condition for the flat model is obtained as 

^2 = (1 - b)H 2 (48) 
Using Eq. @, Eqs. (g) and @ reduce to 

8tt(p - £u*i) = <j)bH 2 (49) 

87rp = -ZcbbH 2 , where b < (50) 

In this case the model has p < which might describe the very early epoch of galaxy 

formation from the process of matter condensation. 

Case (i): b ^ -1. Using Eq. © in Eqs. © and © yield 



^{p-— TTTT^) — 9 (51) 

(l + b)H 2 > {l + b) 2 t 2 

8np = 3 ^ 9 , with b < (52) 

(1 + 6) 2 t 2 



where is already given by Eq. (g9J). Here, the model describes the early phase of evolution 
as mentioned earlier. 
Modell(£ = £o) 



In this case, we assume n = 0. Further from Eq. (34), we obtain £ = £o= constant and 
hence Eq. (BIT) reduces to 



^=(TT^[M^ + 3e ° ] (53) 



Model II (£ = £ p) 



In this case, we assume n = 1. Further from Eq. (|34|), we obtain i = £,oP an d hence Eq.(pT|) 
reduces to 

= ^ b n - i r ™ 

P 8vr(l + 6) 2 i 2[ (1 + &)t J P ; 



Case (ii): = —1. In this case using Eq. (E4|) in Eqs. (49) and (50) yield 



A. Pradhan and H. R. Pandey 9 



8n(p — 3to2^) = — (55) 



8-7TP = 307712 (56) 



where, the Barber's scalar function is given by the Eq. (|39|). 
Model I (£ = Co) 



In this case, we assume n = 0. Further from Eq. (34), we obtain £ = £0= constant and 
hence Eq. reduces to 

p=-^ + 3m 2 Co (57) 

Model II (£ = Cop) 

In this case, we assume n = 1. Further from Eq. (|34|), we obtain £, = CoP an d hence Eq. 
(p5b becomes 

■ — ■ / 2 

p = _^2 (l _ 9m ^ o) (58) 

4. Conclusion 



In this paper we have obtained a class of LRS Bianchi Type I models in Barber's second 
self creation theory in the presence of bulk viscous fluid for constant deceleration parame- 
ter. Assuming an ad-hoc law for the coefficient of bulk viscosity of the form = £, p n , 
where p is the energy density and n is the power index, we have obtained exact solutions of 
the field equations. The nature of the Barber's scalar function <f> and energy density p have 
been examined for values of n = and n = 1 for both the (i) power-law and (ii) exponential 
expansion of both the non-flat and flat universe. The models discussed here are isotropic 
and homogeneous and, in view of the assumption of isotropy, the shear viscosity is absent. 
The effect of bulk viscosity is to produce a change in the perfect fluid. We observe here 
that Murphy's conclusion about the absence of a big bang type singularity in the finite past 
in models with bulk viscous fluid is, in general, not true. 
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